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NORMALIZATION  OF  SPECTRUM  ^ANALYZER  JpUTPUT  DATA  #1 


INTRODUCTION 


If. 


In  general,  the  broadband  continuum  noise  background 
consists  of  either  self  noise,  shipping  noise,  or  ambient  noise  or 
combinations  thereof.  In  any  case  the  spectrum  of  this  noise  will 
not  be  flat,  moreover  the  slope  of  the  spectrum  may  be  highly 
variable  depending  on  the  source  of  the  noise  and  the  range  of 
frequencies  being  considered.  This  condition  is  undesirable  from 
a broadband  detection  point  of  view,  but  from  a narrowband 
detection  standpoint  it  is  an  intolerable  s it.ua tion__if...lef t as  is. 
That  this  is  so  can  be.  seen  readily  by  observing  the  spectral 
representation  in  Fig.  1. 
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Fig.  1 REPRESENTATION  OF  A TYPICAL  BACKGROUND  f.  -/ 
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This  figure  depicts  a situation  where  a narrowband  signal  with  a 
considerably  large  signal- to-noise  ratio  would,  with  high  prob- 
ability, go  undetected  because  of  a non  uniform  power  distribution 
in  frequency.  If  an  OR-Gate  were  employed,  this  loss  in  detect- 
ability would  result  from  capture  of  the  OR-Gate  by  the  frequency 
analyzer  channels  in  the  neighborhood  of  the  frequency  that 
coincides  to  the  maximum  in  the  power  spectral  density  function. 

This  situation  can  be  rectified  by  normalizing  the  spectrum 
analyzer  channel  waveforms  either  before  or  after  envelope  detection 
has  occurred.  Since  we  are  constrained  to  work  with  a hybrid 
spectrum  analyzer  which  does  not  perform  normalization  prior  to 
detection,  it  is  necessary  to  normalize  the  data  after  detection. 

The.  means  by  which  this  is  being  accomplished  is  described  in  the 
remainder  of  this  note. 

PAST  DETECTION  NORMALIZATION 

One  of  the  primary  objectives  of  any  normalization 
scheme  is  to  make  the  probability  of  exceeding  the  threshold 
independent  of  frequency.  This  may  be  referred  to  as  holding  the 
single  channel  false  alarm  probability  constant.  If  we  assume  in 
the  presence  of  noise  alone  (excluding  self  noise  lines)  that  the 
spectrum  analyzer  channel  outputs  are  Rayleigh  distributed,  then 
the  single  channel  false  alarm  probability  is  given  by 


PFA'  J 


2 a 


n 


dx 


(1) 


n 


Now  it  is  clear  that  this  probability  depends  on  only  one  parameter, 

namely  an  . However,  this  parameter  is  the  noise  power  in  the 

bandwidth  of  the  analysis  filter  and  hence  may  be  a function  of 

frequency.  There  is  a way  around  this,  however,  since  the  mean 

2 

standard  deviation,  and  mean  square  of  x are  related  to  by 
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Mx  ■ °nV2 


n 


(2) 


°x  = anV2"-  ^2 


(3) 


E[x2]  = 2c  2 


n 


Now  suppose  that  we  estimate  E[xz]  = P , by  means  of  the  following 
» X 

statistic , 


M 


1 n t 

p = ^ y x.' 

x M l 


and  divide  each  sample  of  x by  P , then  since 

X 


E[P  ] = P 
x x ’ 


we  would  have  a new  random  variable,  y,  given  by 


(4) 


which  for  large  M would  have  a Rayleigh  density  function  with  a 

mean  square  value  of  unity.  If  the  mean  square  of  y is  unity, 

then  the  parameter  in  the  distribution  function  of  y would  be  1/2 
o 

instead  of  a , thus 
n ’ 

' r“  -v2 

pfa  ~ J 2y  e y dy  • (5) 
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2 

This  equation  is  independent . of  the  parameter  on  and  is  thus  not 
dependent  on  the  noise  power  at  the  input  to  the  detector  --  hence 
the  system  is  normalized. 

Now  let  us  examine  exactly  how  this  normalization  is 
being  accomplished  in  the  computer.  The  output  of  the  spectrum 
analyzer  is  digitized  and  entered  into  the  computer  via  digital 
tape  in  a form  which  is  most  conveniently  represented  as  a matrix 
x . 

In  the  time  interval  [t^,  t^]  x has  a form 


x = 


X11  x12 


21 


x 


1L 


m 


NL 


where  each  of  the  correspond  to  the  output  of  the  j frequency 

channel  at  the  t^  instant  of  time. 


that 


To  normalize  the  matrix  x we  define  an  operation  such 


^Lj  ' VPLj 


where 


PT  • “ 
Lj 


i=L-Si 


j-s. 


j+S.+J 


x2  ' 

Xij 


2 , 


+ / XT 

'—i  1J 


j j '=j+s. 
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th 

number  of  frequency  samples  adjacent  to  the  j 
sample  that  are  skipped 


J = number  of  samples  that  are  used  in  one  frequency 
window 


= number  of  time  samples  that  are  used  in  the  time 
window . 


Presently,  a total  of  90_samples  are  being  used  to  compute  Pjj  so 
that  it  is  possible  to  get  normalization  of  the  background  noise  to 
within  i 1 dB  after  transient  effects  have  died  away. 

The  width  of  the  normalization  windows  inthe  frequency 
domain  is  presently  5 Hz  on  either  side  of  the  sample  being 
normalized.  This  means  that  only  non- stationary. ties  that  are  wider 
than  this  will  be  normalized.  Similarly,  in  the  time  domain  the 
normalizing  windows  are  9 seconds  in  duration.  As  a result  we 
expect  to  get  temporal  normalization  only  for  those  waveforms  that 
exhibit  non- stationarities  greater  than  approximately  9 seconds  in 
duration . 


The  normalizer  has  been  checked  out  using  an  array  of 
numbers  with  a non-zero  slope  in  both  the  time  a frequency  domains 
and  it  was  found  to  produce  an  array  of  numbers  all  of  which  are 
constant  except  for  those  around  the  edges  of  the  matrix.  Presently, 
we  are  preparing  to  check  the  normalizer  out  using  sea  data. 
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TKACOR  Project  002  135  01  . 12  December  1969 

ANALYSIS  OF  THE  SIMPLE  OR- GATING  DATA  REDUCTION  TECHNIQUE 
Introduc tion 

It  is  worthwhile  to  develop  theoretical  expressions  for 
the  false  alarm  and  detection  probabilities  of  the  data  reduction 
technique  known  as  OR-Gating  or  the  ^maxdTmun-nf 11  rri  t-pri no  With 
the  expressions  in  hand,  it  will  be  possible  to  predict  the 
performance  of  this  processing  scheme  so  that  comparisons  can  be 
made  not  only  with  experimental  results  obtained  using  the  OR-Gate, 
but  also  with  the  performance  of  other  data  reduction  techniques. 

It  is  the  purpose  of  this  paper  to  develop  such  theoretical 
expressions . 

Analysis 

The  function  of  the  OR-Gate  is  to  examine  N (virtually) 
simultaneous  spectrum  analyzer  output  samples  and  then  to  select 
(gate)  the  maximum.  Symbolically,  if  we  denote  the  OR-Gate  output 
at  time  t^  by  y(t^)  = y j > we  have 

yi  = max  {x.^3,  i = 1,  2,  ...,  N.  (1) 

In  order  to  predict  false  alarm  and  detection  probabilities  we 
must  compute  the  probability  that  y exceeds  a threshold  T when 
noise  alone  is  present  and  when  signal-plus-noise  is  present.  To 
do  this  we  proceed  as  follows.  The  distribution  function  of 
y,  Fy  00  > is  given  by 


u 

l 


1 


F (Y)  = PrCy  < Y], 


but 

y = max  fx^},  i = 1,  2,  ...N, 

so  that 

Fy(Y)  = Pr[max{xi}  < Y] . 

Now  if  the  maximum  of  the  set  {x. } is  less  than  Y,  it  follows  that 
all  elements  of  the  set  are  less  than  Y,  hence 

Fy(Y)  = Pr[xx  < Y,  x2  < Y,  • • ■ XN  < Y] . 

If  the  x^  are  independent,  which  we  assume  here,  then 

Fy(Y)  = Pr[x^  < Y]  • Pr[x2  < Y]  Pr[xN  < Y^> 

but  by  definition  the  distribution  function  of  x,  F (X)  is  given  by 

X 

Fx(X)  = Pr[x  < X]. 

Thus , 

N 

Fv00  “ Fx  (y). 
y i=l  xi 

If  all  of  the  x^  are  identically  distributed,  which  they  will  be 
in  the  absence  of  signal, and  if  the  continuum  noise  is  uniformly 
distributed  (same  distribution  form)  in  frequency,  and  if  the 
normalizer  performs  properly,  then 
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F (Y)  = F V)  . 

yv/n  xv/n 


(2) 


If  one  of  the  channels  contains  signal  plus  noise,  then  F (Y) 
becomes , 


Fy<Y>s+n  ' Fx<YW 

The  probability  of  false  alarm  is  given  by, 


(3) 


PFA  * 1 - VT>n 


- 1 - FxN<T>n  ■ 


while  the  detection  probability  is  given  by 


(4) 


PD  " 1 - VVn 


‘ 1 - Fx(T)s+n. 


(5) 


If  we  define  an  elementary  (single  channel)  probability  of  false 
alarm  by 


PFA  - 1 - Fx<T>n> 


(6) 


and  similarly  a single  channel  probability  of  detection  by 


PD  " 1 ' Fx<TW 


(7) 


then  by  substituting  Eqs.  (4)  and  (5)  we  get 
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- 1 - a - Pp.jN, 


pD  - 1 - a - i’fa>n"1  a - V 


These  equations  reveal  that  the  detection  and  false  alarm 
probabilities  at  -the  output  of  the  OR-Gate  can  be  expressed  in 
terms  of  el  pmcn  taryLq>robaj>illt^a&-.ai:-.the--inpjat---l:n- -the  OR-  Gate  . 
Specifically,  the  probability  of  exceeding  the  threshold,  T,  at 
the  output  of  the  OR-Gate  is  given  by  the  probability  that  at 
least  one  of  the  channels  exceeds  the  same  threshold  at  the  input 
to  the  OR-Gate.  Thus,  if  we  can  specify  the  density  or  distribu- 
tion at  the  output  of  the  analyzer  channels,  we  can  get  the 
required  probabilities  of  false  alarm  and  detection. 

When  noise  alone  is  present  and  when  an  envelope 
detector  is  used,  we  have  as  the  density  function  of  x, 

X2 

‘ 2a  ? 

p(x)  = — e , 


where 


on  = noise  power  in  each  analyzer  channel 


Thus , 


J —2  e 

T °n 


which  is 


. -T2/ 2ci  2 
PFA  " e 
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* 

When  signal  plus  noise  is  present  in  the  channel,  we  get  for  P^, 

. x2 

, » 2 (a  2 + o 2) 

PD  = ^ 2~  2 e S n d- 

T s n 

This  assumes  that  the  signal  is  a narrowband  Gaussian  process  with 
2* 

power  o's 

Letting, 


gives 


PD  I 1+a 


e dy 


where 


a.=  , the  signal- to-noise  ratio. 

a 1 
n 

Equation  (11)  can  be  integrated  by  letting  w = j thus 


Pn  = e 


2o/(l+a) 


This  is  a reasonable  assumption  when  one  considers  the  effects  of 
the  ocean  channel.  That  is,  a pure  tone  originating  at  a source 
will  undoubtly  suffer  considerable  fading  due  to  multipath  effects 
and  will  result  in  a received  signal  that  is  the  superposition  of 
many  "tones"  which  have  random  amplitudes  and  phases  with  correla- 
tion time  comparable  to  reciprocal  bandwidth. 
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65C 

By  defining  a normalized  threshold, 


T = 


2 a 


n 


Equations  (10)  and  (12)  become 


-T 


PFA  = C 


P'  = -T'/d+a) 

rD  6 


(13) 


(14) 


Using  Eqs.  (13)  and  (14)  we  can  get  the  required  false  alarm  and 
detection  probabilities,  thus 


rp  i'l 

Fa  = 1 - (!  “ e ) 

(15) 

-T')N'X  (1  - e-T'/<1+a>). 

(16) 

D 


With  Eqs.  (15)  and  (16)  it  is  possible  to- -compute- .ROC  curves  for 
the  simple  OR-Gate  as  _a  . function  of  the  s-i  gna1  - tn-noise  ratio,  a. 
These  curves  may  be  used  for  comparison  with  experimental  results 
as  well  as  for  comparing  different  data  reduction  techniques. 


ANALYSIS  OF  THE  MAXIMUM  LIKELIHOOD 
RATIO  DATA  REDUCTION  CRITERION 

INTRODUCTION 

The  maximum  likelihood  technique  for  deciding  which 
spectrum  analyzer  channel  is  most  Likely. to  contain  a target 
spectral  line  is  based  on  computing  the  likelihood  ratio  of  the 
data  at  the  output  of  each  channel  of  the  spectrum  analyzer  and 
then  subjecting  the  set  of  all  such  likelihood  ratios  to  an  0R- 
Gate.  That  is,  the  output,  y(t),  of  such  a system  is  given  by 

y(t)  = max  {L^x^t)]}  , (1) 

i 


where 


i = 1 , 2,  . ..,  N = the  spectrum  analyzer  channel 
number,  and 

th 

x^(t)  = the  time  function  at  the  output  of  the  i 
spectrum  analyzer  channel. 


P(x-  )s-ln 
L.  (x.  ) = ---/  v- 

1 1 p^xiV 


(2) 


where 

Ps+n  (xd  ) = Pr°bability  density  function  of  the  it^1 

spectrum  analyzer  channel  when  signal  plus 
noise  is  present  at  its  output,  and 

p(x^)  = probability  density  function  of  the  i*7*1 
spectrum  analyzer  channel  output. 

There  are  two  interpretations  that  can  be  placed  on  the  likelihood 
ratio,  the  first  of  which  is  that  it  is  the  likelihood  ratio  for 
single  output  sample;'..  The  second  interpretation  is  that  L(x) 
is  a joint  likelihood  ratio  obtained  by  operating  on  successive 


1 


t 


\ 

i 
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outputs  of  each  analyzer  channel.  It  is  the  purpose  of  this  note 
to  examine  these  possible  implementations  and  the  consequences  of 
each. 


SINGLE  SAMPLE  LIKELIHOOD  RATIO 


In  order  to  establish  this  likelihood  ratio  it  is  neces- 
sary to  first  determine  the  form  of  the  two  required  density  func- 
tions, p(x)g+n  and  p(x)n.  The  function  p(>0n  will  be  a Rayleigh 
density  function  if  the  spectrum  analyzer  output  detector  is  an 
envelope  detector  and  if  the  analysis  filter  output  is  Gaussian. 
Both  of  these  conditions  are  reasonable;  the  first  because  of 
practical  considerations  and  the  second  because  of  physical  con- 
siderations. The  Rayleigh  density  function  is  given  by 


P 


e 


(3) 


where 

2 

= the  noise  power  level  at  the  input  to  the  detector, 
i.e.,  in  the  bandwidth  of  the  analysis  filter. 

To  obtain  the  density  of  the  detector  output  in  the  presence  of 
signal  plus  noise  it  is  necessary  to  first  define  what  is  meant 
by  signal.  The  most  reasonable  model  for  a narrowband  signal  is 
that  of  a narrowband  Gaussian  process.  This  is  true  because  of 
the  multipath  effects  that  will  be  experienced  by  tone- like  sig- 
nals in  traveling  from  the  source  to  the  receiver.  Thus  the  input 

to  the  envelope  detector  can  be  modeled  by  the  sum  of  two  Gaussian 

2 

processes,  one  (the  signal)  with  power  a and  the  other  (the  noise) 
2 . s 

with  power  . This  means  that  the  envelope  detector  output  is 
described  by  another  Rayleigh  density  function,  i.e., 


P(>0 


s+n 


2/o/  2,  2. 

-*  /2("s  +0n  > 


O) 
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The  likelihood  ratio  is  given  by  the  ratio  of  Eqs.  (4) 
and  (3) , thus 


L(x)  = — j 


n 


+(x2/2an2)-[x2/2(as2-fon2)] 
e , (5) 


a + a 
s n 


Since  the  logarithm  is  a monotonic  function,  we  can  just  as  easily 
consider  the  logarithm  of  the  likelihood  ratio,  i.e., 


t(x)  ^ log  L(x) 


where 


C1  = 


% (1  + a) 


= C-.X  + C9  , 


a 

~~7~2 7 ’ 

Gs  + an 


(6) 


Co  = logo(- 


e.  a + 1' 


and 


a = 


2 

°s 
— 7 
°n 


the  pre- detector  signal- to-noise  ratio,  this 
is  a so-called  design  signal- to- noise  ratio. 


Some  interesting  observations  can  be  made  concerning 

Eq.  (6).  First,  in  order  to  compute  -t(x)  we  have  to  know  or 
. 2 

estimate  both  the  noise  power  an  and  the  signal- to-noise  ratio, 
a.  This  means  that  we  would  have  to  design  for  some  minimum  value 
of  "a"  and  estimate  on  , the  latter  being  a simple  job  if  the 
normalization  is  done  properly.  Second,  it  is  clear  that  -t(x) 
is  a monotonic  function  of  x and  thus  if  the  i^1  channel  output, 
x^,  is  the  maximum  of  all  the  frequency  channels,  then  it  follows 
that  L^(x^)  will  also  be  the  maximum  over  all  the  likelihood  ratios 
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of  the  channel  outputs.  That  is  to  say,  if  the  i1"'1  channel  is 
chosen  by  a simple  OR- Gating  of  the  frequency  analyzer  channel 
outputs,  then  the  iu  channel  will  also  exhibit  the  maximum  like- 
lihood ratio  and  will  thus  be  chosen  again.  The  next  question  to 
be  considered  is  whether  this  maximum  likelihood  ratio  technique 
will  give  the  same  probability  of  detection  for  equal  probability 
of  false  alarm.  To  investigate  this  we  proceed  as  follows. 

Let  the  output  of  each  frequency  analyzer  channel  be 
distributed  according  to  Eq.  (3)  in  the  presence  of  noise  alone 
and  according  to  Eq.  (4)  in  the  presence  of  signal  plus  noise. 

Let  the  false  alarm  probability  be  given  by  the  probability  that 
the  output  of  the  OR- Gate  exceeds  a threshold  T.  That  is, 


PFA  " 1 - Fy<T>  • 


(7) 


where 


Fy(T)  = Pr(y  < T) 

y = max  [xi (t) ] . 

The  distribution  function  of  y is  given  by 


Fy (T)  = FxN(T)  , 


(8) 


where 


FX(T) 


We  note  that 


- Pr(x  < T)  = distribution  function  of  each  of  the 
identically  distributed  x^'s. 


Fx(T)  = 1 - Pr(x  > T) 


= 1 _ p 7 

L rFA  » 


(9) 
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fnr-wmrri 

where 


Pp^  = probability  of  false  alarm  in  each  spectrum 


analyzer  channel. 

Substituting  Eqs.  (9)  and  (8)  into  Eq.  (7),  we  get 


PFA  1 ~ (1  " PpA) 


N 


(10) 


This  equation  measures  the  probability  that  at  least  one  of  the 
channel  outputs  exceeds  the  threshold  T.  What  has  been  shown  is 
that  we  can  predict  performance  at  the  output  of  the  OR- Gate  by 
performing  the  analysis  at  the  input  to  the  OR- Gate.  This  approach 
agrees  with  our  intuitive  reasoning  if  we  interpret  the  problem  as 
one  of  thresholding  before  OR-Gating. 

Substituting  Eq.  (3)  into  Eq.  (10)  we  get. 


PFA  1 


a - J' 


2/9  2 

-x  /2a  v, 

x e n dx)N 

°n 


(ID 


for  an  explicit  form  of  the  false  alarm  probability. 

The  integral  in  Eq.  (11)  can  be  evaluated  easily  and  is 


given  by 


! 


" X -*2/2on2  -T2/2o  2 

T —j  e dx  = e 

an 


(12) 
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The  detection  probability,  P^,  is  given  by  the  probability 
that  when  a signal  is  present  in  one  of  the  channels,  that  the  0R- 
Gate  output  will  exceed  the  threshold  T.  That  is, 


where 


PD  = 1 - (1  - PfA)N_1(1  - Pp)  , (14) 


-x2/[2(as2+an2)] 


dx  . 


(15) 


But  from  the  preceeding  agreements  we  have 

-T2/[  2 (o  2-fcr  2 )] 
P ' = e sn 


(16) 


Substituting  Eqs.  (16)  and  (12)  into  Eq.  (15)  gives, 


(1  - e 


•T2/2o 


n xN-1 


)1N-X(1 


- T2/[ 2 (ct  2-fo  2 )] 
e S n ) . 


(17) 


We  have  thus  expressed  the  false  alarm  probability  and  detection 
probability  for  the  case  where  the  spectrum  analyzer  output  is 
directly  driving  the  OR- Gate.  What  must  be  done  next  is  to  show 
that  when  the  single  sample  likelihood  ratio  is  employed,  the 
probability  of  detection  is  the  same  as  that  given  by  Eq.  (17) 
when  the  false  alarm  probability  is  the  same. 


If  in  forming  the  logarithm  of  the  likelihood  ratio  x 

is  mapped  according  to  Eq.  (6),  then  the  probability  density  func- 

2 

tion  of  the  new  random  variable,  z = -t(x)  = C-^x  + C2 , is  given  by 


P(z) 


with  x 


1 
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Since 


2Clxdx> 


2 - Cr 


" C 2 ’ 


and  given 


2/0  2 
-x  2 a 

/ \ x ' n 

p(x)  = j e 


■we  get  for  p(2), 


l -z-C2/2C1Gn" 

p(z)  = 7 e 

2clCTn 


The  probability  of  false  alarm,  Pp^p  for  this  case  is 


FFA1  - 1 - <!  - j' 


1 -z-C2/2c  ]_an^ 


T1  20^ 


dz)  , 


or  in  simplified  form 


PFA1  - 1 - (1  - eC2/2Cl""2  . 
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By  similar  reasoning  we  can  .show  that  the  detection  probability 
PD1  is  given  by 


PD1  " 1 " (1  - P.ai) 


N-l 


FA1' 


a - J' 


T1  2Cl(a/  + an‘) 


1 -z-C2/[2C1(as2+an2)]  ^ 

77  e 


or 


+C9/2C1an2  -T,/2C,an2  „ , 

PD2  = 1 - (1  - e 2 1 n e 1 1 n )N"1 


(19) 


(1  . e+C2/C2c1<as2+CTn2)3  e-T1/[2C1(os2+an2)]  • 


Equations  (17)  and  (19)  give  the  required  probabilities 
of  detection  as  functions  of  the  probabilities  of  false  alarm,  so 
that  what  remains  to  be  shown  is  that  for  equal  false  alarm  prob- 
abilities we  get  equal  detection  probabilities. 

If  PFA  = pFAi>  this  means  that,  using  Eqs.  (13)  and  (18) 

T2/2a  2 . . 

2C-.0  2 
1 n 


or  that  in  order  for  our  assertion  to  hold  it  must  be  true  that 


8 


which  reduces  to 


(22) 


Substituting  the  right  hand  side  of  Eq.  (22)  for  in  Eq.  (21) 
and  simplifying  gives 


which  proves  the  assertion. 

What  we  have  demonstrated  is  the  rather  intuitively 
obvious  fact  that  by  performing  a monotonic  mapping  on  the  out- 
put of  the  spectrum  analyzer  channels,  the  probability  of  detec-; 
tion  for  a given  false  alarm  probability  is  not  changed.  Thus 
it  may  be  concluded  that  there  is  nothing  to  be  gained  by  per- 
forming  such  an  operation. 
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JOINT  LIKELIHOOD  RATIO 

If  wc  make  the  same  assumptions  concerning  the  target 
signal  and  background  noise,  i.e.,  that  they  are  Gaussian  processes, 
then  it  is  possible  to  compute  the  joint  log  likelihood  ratio  for 
the  output  of  each  analyzer  channel.  The  joint  likelihood  ratio 

A.V. 

for  M independent  1 analyzer  channel  outputs  is  given  by 


where 


. , , e<xlVn 
l&i>  - 57577“  > 


[x^ (^"i^  » xj_  ’ • • • > 3 


P s+n  P^xil^x+n  P^xi2^s+n  ""  P^xM^s+n 


P(x)  = P(xil)np(si2)n  •••  p(xiM)n 


(23) 


That  is,  Eq.  (23)  can  be  expressed  as 


P(xn) 

s+n  . P<xi2>s+n 


T / \ J-JL  o in 

<£i>  ' P(xil>n  'P^‘i2>n 


P^xlMh4n 

' P^Wn 


If  we  take  the  logarithm  of  the  above  equation  and  use  the  earlier 
l(x)  = logeL(x),  we  get 


M 


UiSi)  = Y,  ^ji)  * 

j=l 


(24) 


10 
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Recalling  from  Eq..(6)  that 


-t(x)  ~ Chx^  + C0  , 


we  get  for  Eq.  (24) 


4(2^) 


M 

Cl  I x 

3=1 


+ NC, 


(25) 


This  equation  tells  us  that  the  process  under  consideration  is 
square  law  detection  and  summation  of  M variates.  The  next  step 
is  to  determine  the  probability  density  functions  of  4(x^)  under 
conditions  of  noise  alone  and  signal  plus  noise.  Define  y as 


M 


* 


a - * E xfi  • 


- 2 

Of  X.. 


When  signal  plus  noise  is  present  we  have  for  the  density 


--  2 


p(z)  = _1 2 


-z/(CTs2+a n2> 


(26) 


CT  + 

s n 


The  factor  \ is  used  here  for  mathematical  convenience  and  is 
equivalent  to  using  a square  law  detector  with  a gain  of 


11 
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.1,2 


Others  ’ have  shown  that  the  density  of  is  given  by 

it  2.  2. 

i yn-  m i -y< /(<*«:  +CTn  ) 

~ (M  - 1) ! 7— — 2>  e • (27) 

as  an 

2 

To  get  the  density  function  for  noise  alone  we  let  ag  =0  and 
Eq.  (27)  becomes 


, . 1 ^isM-l  -?i/an 

TM^-I1T(7T’  e 


(28) 


n 


The  false  alarm  probability  after  OR-Gating  is  given  by  Eq.  (10) 
and  when  Eq.  (28)  is  substituted  we  get 


PFA  ~ 1 " C1  -J  T IM  ^LnT("J12) 


M-l 


~y/on  dy  vN 


-^7).  (29) 


— 7. 

cn 


n 


n 


Similarly,  the  detection  probability  is  given  bv 


PD  = 1 " (1  "I  T '(M'"-1'T)T(f7)M’1  e n -^-)N'1 


n 


n 


n 


(30) 


" -r. 


1 -y/(°s2+Pn2) 


T 


T 

""  "2.  2 

CT 

s n 


a + a 

s n 


27  2 
a +CT 
s n 


Marcum,  J.  1. , "A  Statistical  Theory  of  Target  Detection  by 
Pulsed  Radar:  Mathematical  Appendix,"  RM-753,  1 July  1948, 

The  Rand  Corporation. 
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Equations  (29)  and  (30)  can  \>e  expressed  in  terras  of  a normalized 
threshold  thus  giving 


FA 


- 1 


- a - j; 


1 M- 1 -w  , 

IK  1') ! w e d“> 


N 


(31) 


and 


where 


u — i ( i r 1 M—  1 - w j \ N- 1 

pd  " 1 ■ (]  - Jt7  ~(M~iyr  w e dw) 


(i  - 


r i , sM-i  -w  , .. 

j , ~(m  - TJT  'w)  e dw) 

_T_ 

1+a 


T7  = T/u  1 , and 


(32) 


SUMMARY 


g 

a = — y , the  signal- to- noise  ratio. 
an 


a s 


We  have  developed  two  results  lay_tlie  analysis,  developed 
in  this  paper.  First,  it  s ho wn  that  nothing  is  to  be 

gained  by  single- sample  likelihood  ratio  processing.  Although 
this  is  an  intuitively  apparent  fact  it  never  hurts  to  verify  such 
assertions.  The  other  result  is  that  it  is  posslbleto  tlieoreti- 
cally  predict  the  performance  of  the  joint  likelihood  ratio,  data 
reduction,  criterion  under  quite  realistic  assumptions  concerning 


4^1-)  and  (32)  can  be  cval- 


the  noise  and  signals.  There fore^ 
uated  numerically  and  receiver  operating  characteristic  (ROC) 
curves  can  be  obtained.  These  curves  can  be  compared  to  the  (ROC) 
curves  that  are  obtained  both  theoretically  and  experimentally 
Trom  the  simple  OR- Gating  data  reduction  technique. 


